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The quarter-filling extended Hubbard model on the triangular lattice is studied to explore 
pairing instability in the three-fold charge-ordered (CO) metal. We derive a second-order strong- 
coupling effective Hamiltonian of doped carriers into the three-fold CO insulator at electron 
density of n = 2/3, and then study the /- and d^-wave superconductivities down to n = 1/2 
by using the BCS mean-field approximation. It is found that the triplet /-wave pairing is more 
stable than the d^-wave one. We also discuss that this coexisting state of the charge ordering 
and superconductivity is possible to have critical temperature T c ~ O.Oli. 

KEYWORDS: triangular lattice, three-fold charge order, extended Hubbard model, superconductivity, BCS 
mean-field approximation, organic conductor, #-(BEDT-TTF)2X 



Motivated by experimental investigations on (9-type 
organic conductors 1 and Na^Co02-yH 2 0, 2 the possi- 
ble relevance of charge ordering and superconductiv- 
ity has been studied extensively on the basis of the 
triangular-lattice extended Hubbard model with an ad- 
ditional nearest-neighbor Coulomb repulsion V. 3 ~ 7 

Charge-fluctuation mediated superconductivity in the 
vicinity to a CDW phase is one of main topics in those 
investigations. Tanaka et al. applied the RPA method to 
the triangular-lattice extended Hubbard model, having 
the superconductivity in Na x Co02-yH20 in their mind. 3 
They found that V enhances the charge fluctuation at 
q = (2tt/3, 27r/3) and triplet next-nearest-neighbor /- 
wave superconductivity is induced. Onari et al. obtained 
phase diagram by using the FLEX approximation. 4 They 
found that /-wave pairing adjacent to a CDW phase is 
stable for intermediate filling and large V. These results 
have been paid attentions, because triplet superconduc- 
tivity was considered to be stabilized only under some 
special conditions. 

Due to frustration effects inherent to the triangular- 
lattice, a variety of charge-ordering patterns have been 
considered in this system. 5-10 Watanabe et al. stud- 
ied the triangular-lattice extended Hubbard model at 
n = 2/3 by the variational Monte-Carlo (VMC) method, 
and found two kinds of CO states, three-fold bipolaronic 
CO state and antifcrromagnctic CO state, in the U-V 
phase diagram. 5 In the three-fold CO state, A-sublatticc 
sites are doubly occupied and B- and C-sublattices have 
no electrons, as shown in Fig. 1(a). In the antiferromag- 
netic CO state, there exist no electrons on A-sublattice. 
Each of B- and C- sublattice sites is singly occupied, and 
the spins are in honeycomb- type antiferromagnetic ar- 
rangement. Ground-state energies per site of the three- 
fold and antiferromagnetic CO states are, respectively, 
U/3 and V in the strong-coupling limit, and thus we ex- 
pect the transition point V ~ U/3. And then they stud- 
ied an anisotropic-triangular-lattice extended Hubbard 
model at quarter filling, which is a model for #-(BEDT- 
TTF) 2 X, by the VMC method. 6 They found that there 
appears three-fold CO metallic phase rather than insu- 



lating stripe CO phases if the anisotropy in the intersite 
repulsions is not so large. 

Organic conductors 6»-(BEDT-TTF) 2 X have been 
known to show superconductivity when X=l3. Watan- 
abe et al. speculated that 6»-(BEDT-TTF) 2 I 3 is located 
in the vicinity of the para metal phase in their ground- 
state phase diagram, 6 and thus it is a weak coupling sys- 
tem. Aside from the description of organic conductors 
by a realistic theoretical model, there arises a fundamen- 
tal question: Is the extended Hubbard model possible to 
show superconductivity in intermediate and strong cou- 
pling regions? When Coulomb interactions are increased, 
then a charge ordering is expected to be more favored. 
If the CO state has a charge excitation gap, then super- 
conductivity is never obtained. On the other hand, if the 
CO state is metallic, then charge carriers are possible to 
be condensed into a superconducting state. Such super- 
conducting state is inevitably coexisting with the charge 
ordering. 

The purpose of this paper is to explore theoretically 
the above-mentioned coexisting state in the triangular- 
lattice extended Hubbard model with quarter-filled 
band. To be precise, we study pairing instability in the 
three-fold CO metallic phase. Unfortunately, it has not 




Fig. 1. Schematic representations of (a) the three- fold CO insu- 
lator realizing for V>[//3atn = 2/3 and (b) the three-fold 
CO metal which is obtained by hole doping into the CO insu- 
lator. Second-order perturbation processes make doped holes be 
itinerant within A-sublattice, which is schematically represented 
by the gray arrows. 
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been reported that the three-fold CO metallic phase is 
realized in organic conductors. We hope that the three- 
fold CO metallic phase is found experimentally in the 
near future, which opens new windows to the physics of 
charge ordering and superconductivity. 

In order to study the present issue, we use the extended 
Hubbard model on the isotropic triangular-lattice, 

H = -t ^ (c^Cja- + h.C.) + U~^2 n i1 n il + V E niTl ^ 

(1) 

where denote nearest- neighbor pairs of lattice 

points. Hereafter, we assume the strong-coupling limit, 
( « [/, V, in our calculations. For the stability of the 
three-fold charge ordering, we also assume U < 3V. 

Our starting point is the three-fold CO insulator at 
n = 2/3, which is schematically shown in Fig. 1(a). If 
we make hole doping to this insulator, the three-fold CO 
metal is obtained. Note that the idea of carrier doping 
to an insulator is similar to the high-T c cuprates. 11 Fig- 
ure 1(b) illustrates the strong-coupling picture of how 
doped holes propagate in the CO background: a hole can 
hop from an A-sublattice site to another A-sublattice site 
next to the first one via the second-order perturbation 
process of the bare hopping t. 

Some readers may wonder if the hole doping into a 
bipolaronic CO insulator collapses the charge order dras- 
tically, like doped antifcrromagncts. 12 This issue has 
been studied by several authors in the linear-chain, 13 
two-leg ladder, 14 and square lattice. 13, 15 Those studies 
indicated that charge order is robust against hole doping 
in the intermediate- and strong-coupling regions. As for 
the triangular lattice treated in this paper, Hotta et al. 
recently studied the t-V model for spinless fermions at 
the density p ~ 1/3. 8 They confirmed that the three-fold 
charge ordering survives even at p > 1/3 and the charge 
excitation gap closes there. Note that charge order can 
be stabilized only by the Coulomb repulsion, although 
both the hopping and Coulomb repulsion are required 
to stabilize antifcrromagnctic order. The difference be- 
tween CO and antifcrromagnctic states may be coming 
from this fact. 

Now we turn to our theoretical analysis. In order to 
identify an ordering stabilized in the doped three-fold CO 
system, we calculate the second-order effective Hamilto- 
nian within the two- hole approximation. 16 The obtained 
effective Hamiltonian is an triangular-lattice extended 
Hubbard model with additional nearest-neighbor inter- 
action and correlated hopping terms: 

H eS = T 1 Y / hfo +V Q J2 »?T n 2l 

+ V 1 J2^>i'+TiY, E n t" h t^ ( 2 ) 
{W) {i,i')t"eT iA , 

where I denotes an A-sublattice site, (£, £') a nearest- 
neighbor pair of A-sublattice sites, a set of A- 
sublattice sites which form regular triangles together 
with the lattice points £ and £', hp e , = 2 CT ( a L a ^' '+ n - c -) 
with a hole annihilation operator ai a = c\ a (hole-particle 
transformation), n a la = a\ a ai a i n>g = The matrix 



elements in H c g are given by 

2t 2 1 

12t 2 / 1 2 4 4 \ , , 

V ° = U - — (-4 + 3^-4-^+5-W' (4) 

Vl = -— {l2 + J^- — + —)> (5) 

fl = v {jh ih) ' (6) 

where w = U/V. 

Let us examine each of the matrix elements. First, T\ 
makes a hole hop from an A-sublattice site to another 
nearest-neighbor A-sublattice site, as described before. 
Secondly, eq. (4) tells us that the on-site interaction Vb 
is repulsive and Vq ~ U. Thirdly, we find the nearest- 
neighbor interaction V\ being attractive in eq. (5), and 
thus it could be a driving force of superconductivity. 
The origin of this attraction is as follows. We suppose 
that there exists a nearest-neighbor pair of doped holes, 
and consider a process in which an electron hops from 
an A-sublattice site to a B-sublattice (or C-sublattice) 
site next to the doped hole pair (see Fig. 2). The energy 
difference between the initial and intermediate states is 
3V — U, which vanishes in the limit U/V — > 3. Such 
small energy difference is never obtained for more distant 
doped hole pair, so the nearest-neighbor interaction be- 
comes attractive. Finally, the matrix element of the cor- 
related hopping, Ti, is positive. It has been known that 
correlated hopping induces superconductivity depending 
on positive and negative of its matrix element. 16, 17 As 
shown later, the positive T± actually stabilizes the /- and 
d X y-w&ve pairings. 

We here derive the dispersion relation of free holes 
by neglecting the interaction terms in i? c ff- Making 
the Fourier transformation a\~ a = N^ 1 ^ 2 J2 e e lk ' re ai a , 
where N is the total number of A-sublattice sites, 
we obtain the dispersion relation e& = 27\ [cos k y + 
2 cos(v / 3fca;/2) cos(fc y /2)], where the unit of length is the 
lattice constant of A-sublattice. The band width is given 
by 9Ti. 




Initial state E=U Intermediate state E = 3V 

Fig. 2. Schematic representation of a second-order perturbation 
process by a bare hopping term, where the black and gray arrows, 
respectively, represent electrons and holes. This process is the 
origin of the nearest-neighbor attraction between two holes. The 
energy of the initial state is U because of the double occupancy, 
and that of the intermediate state is 3V because of the intcrsitc 
repulsion V. (See text in detail.) 



J. Phys. Soc. Jpn. 



Letter 



S. Morohoshi and Y. Fukumoto 3 



We turn to the study of superconductivity on the ba- 
sis of the obtained effective Hamiltonian. Here we con- 
sider two types of pairing, /-wave and d xy -wa.ve shown in 
Fig 3. We define the off-diagonal mean-field parameter as 
A(p) = (a^at+pi) , where p represent nearest-neighbor 
vectors of A-sublattice. The /- and d xy -~w&ve order pa- 
rameters satisfy the symmetry relations A(p) = — A(— p) 
and A(p) = A(— p), respectively. We also denote the ab- 
solute value of A(p) as A. 




CO-AF boundary 





Fig. 3. Schematic representations of next-nearest-ncighbor (a) /- 
wave and (b) d^y-wave symmetries in real space. In fc-space, the 
gap function changes its sign six (four) times for /-wave (d xy - 
wave). The dotted lines represent the original lattice. 



Applying the BCS mean-field approximation to H c s 
and using the Bogoliubov transformation, we ob- 
tain the diagonalized mean-field Hamiltonian H m f = 
Hk,a E ka\ a ot ka + Efcfe ~ E k) + 7A 2 , where a k<7 de- 
notes the annihilation operator of a quasi particle and 
E k = y/? k + A 2 /fc with i k = e k - p, 



fk = 



ky I "\/3k x ky\ r r 

— sm — cos cos — for r-wave 

3 2 I 2 2 I J 



sm ■ 



\^3k x , k v 

sm — 

2 2 



for d X y-w&ve 



and A = 7A. The pairing interaction 7 is defined by 



(7) 



7 = 



f6(4Ti-Vi) 
[4(2fi-Vi) 



for /-wave 
for d xy -wa,ve 



(8) 



Minimizing the free energy with respect to the order pa- 
rameter A, we obtain the following mean-field equation, 



1 



2tt 2 



f 2 

,^tanh 
2E k 



2k B T 



(9) 



The chemical potential /i is determined by the condition 



which turns out to be 



n h 



V3 
2tt 2 



Jo 



f — — tanh 
E k 



E k 
2k B T 



(10) 

In eqs. (9) and (10), we have transformed the integration 
range in the fc-space from the hexagonal first Brillouin 
zone of the triangular-lattice to the rectangular region. 

In Fig. 4, we show the U /V dependence of the pair- 
ing interaction 7, divided by the band width 9T\, for the 
/-wave pairing. We obtain larger 7 when we get closer 
to the CO-AF boundary U/V = 3. The value of 7 ex- 
ceeds the band width for U/V >~ 1.6. We now com- 
ment on the applicability of the BCS mean-field theory 



Fig. 4. Dependence of the pairing interaction 7 on U/V for the 
/-wave pairing. The CO-AF boundary is [U /V) c = 3. 



to our effective Hamiltonian. In the context of BCS- 
BEC crossover, the relation between critical tempera- 
ture T c and the strength of pairing interaction was ex- 
tensively studied on the basis of the three-dimensional 
negative-// Hubbard model. The calculation of T c for 
this model at quarter filling was carried out by using the 
self-consistent T-matrix approximation 18 and the quan- 
tum Monte-Carlo (QMC) simulation. 19 The result of 
the QMC study is that T c takes the maximum value 
T* ~ 0.35t at \U\ = \U*\ ~ 8t (see Fig. 2 in rcf. 19). 
The BCS result agrees with the QMC result only for 
\U\ -C \U*\. On the one hand, if \U\ exceeds the band 
width, then the BEC formula T c = <b X [numerical con- 
stant (~ 3)], where t B = 2t 2 /\U\ is the effective hopping 
amplitude of composite bosons, becomes to be a good 
approximation. Our pairing interaction 7 is so large that 
we can not use the BCS mean-field theory to give quanti- 
tative predictions. We here use it to get only qualitative 
aspects such as relative stability between the singlet and 
triplet pairings, filling dependence of superconductivity, 
and so on. In the estimation of critical temperature, we 
use the BEC formula instead. 

We turns to the solution of the self-consistent equa- 
tions (9) and (10). In Fig. 5, we show the order param- 
eter Ao = A(T = 0) and the critical temperature T c for 
the /- and d^ y -pairings as a function of hole density of A- 
sublatticc, (= 2— 3n). It is found that the introduction 
of holes leads to rapid development of superconductivity, 
and Ao and T c saturate at ~ 0.5. Thus, the quarter 
filling, n = rih = 0.5, is nearly optimal doping rate for su- 
perconductivity. Although the overall behaviors are sim- 
ilar between those two pairings, T c and Ao of /-pairing 
are larger than those of d^-pairing for all U/V, which is 
mainly coming from the fact that the d^-pairing interac- 
tion is smaller than the /-pairing interaction due to the 
vanishing of two of A(p) in the d^^-wave case. In order 
to confirm the stability of the /-wave superconductivity, 
we calculate the temperature dependence of the free en- 
ergies of the /- and d xy -wave states per A-sublattice site. 
The calculated result for U/V = 2.5 and quarter filling 
is shown in Fig. 6, where we find that the free energy of 
the /-wave state is lower than that of the <i xy -wave state 
for the whole temperature range. Thus, we conclude that 
the /-wave superconductivity is stable in the three-fold 
CO metal. 

Finally, we try to estimate T c in the intermediate 
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Fig. 5. Hole-density dependence of (a) zero-temperature order 
parameter Ao and (b) critical temperature T c for the triplet /- 
wave and singlet d xy -wa.ve superconductivities, where the units 
of order parameter Aq and critical temperature T* are the values 
for the /-wave pairing at quarter-filling, nj, = n = 0.5. 





n = 1/2 ; 




U/V=2.5 i 
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r d -wave 
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Fig. 6. Mean-field free energy F for the /- and d xy -wave pairings 
as a function of T at quarter filling, where the temperature is 
normalized by T c of the /-wave pairing and F* represents the 
value of free energy for the /-wave pairing at T = T c . The dotted 
line represents free energy in case of no pair condensation. 



coupling region, assuming our strong-coupling expansion 
holds even there. As mentioned previously, the estima- 
tion of T c itself within the BCS theory, which gives un- 
reasonable large values, is meaningless. Instead, we cal- 
culate effective hopping amplitude of composite bosons, 
f B = 2^/7, which gives a measure of T c for our effective 
Hamiltonian. Using 7 for the /-wave pairing, we show te 
as a function of on-site repulsion U in the original Hamil- 
tonian in Fig. 7. This result suggests T c ~ O.Oli for the 
superconductivity in the three-fold CO metallic phase. 
Also, our strong pairing interaction may cause phase sep- 
aration near the CO-AF boundary. This speculation is 
based on the study of the extended Hubbard ladder by 
the present authors, 20 where a similar effective Hamilto- 
nian was analyzed by using the Luttingcr theory. They 
found negative Luttinger parameter [K p < 0), which in- 
dicates phase separation, near the CO-AF boundary. 

In summary, we have studied the triangular-lattice ex- 
tended Hubbard model by using strong-coupling expan- 
sion. It has been found that charge carriers in the three- 
fold CO metal shows the triplet next-nearest-neighbor 
/-wave superconductivity with T c ~ O.Olt. On the one 
hand, our theory is based on the strong-coupling effec- 
tive Hamiltonian within the two-hole approximation. It 
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Fig. 7. Effective hopping amplitude of composite bosons, *b = 
27^/7, for the triplet / -wave superconductivity as a function of 
on-site Coulomb interaction U. 



is desired to study how the omitted many-hole processes 
affect the present results. However, the straightforward 
extension of our perturbation theory is difficult and hope- 
less. An alternating way of treating this problem is the 
VMC method. Giamarchi et al. proposed a variational 
wave function with both the superconducting and mag- 
netic orders in their VMC study of the Hubbard and 
t-J models. 21 In a similar way, we can construct a wave 
function with both the superconducting and three-fold 
charge orders. Such investigations are now in progress. 
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